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ABSTRACT 
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appoximating  two-point  boundary  value  problems  for  linear  ordinary 
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SIGNIFICANCE  AND  EXPLANATION 

m 

The  normal  procedure  in  solving  a  continuously  defined  problem 
nwerically  consists  in  applying  a  discretisation  first  which  reduces  the 
original  problem  to  one  which  can  be  treated  by  numerical  algorithms  for 
solving  equation  with  a  finite  number  of  unknowns.  A  fundamental  question 
arising  in  this  context  is  the  convergence  of  the  solutions  of  the  dlscretised 
problem  to  the  eolution  of  the  original  one. 

In  studying  these  kinds  of  questions,  it  has  turned  out  that  the  various 
different  discretisation  procedures  can  be  treated  in  a  unified  manner  in  the 
abstract  setting  of  the  so-called  discrete  approximation  theory  which  is  also 
used  in  this  paper.  The  main  point  in  proving  convergence  in  the  stability  of 
the  discretisation.  This  paper  deals  with  methods  in  proving  stability  for 
linear  problems.  As  an  application  of  the  abstract  result  obtained,  the 
stability  of  finite  difference  aproximations  for  linear  two-point  boundary 
value  problem  on  not  equally  spaced  grids  is  established,  a  topic  which  has 
attracted  the  attention  of  a  number  of  numerical  analysts  in  the  last  few 
years. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
■immirj  lies  with  MAC,  and  not  with  the  author  of  this  report. 


OH  PERTURBATIONS  OF  STABILITY  INEQUALITIES  KITH  AN  APPLICATION 
TO  FINITE  DIFFERENCE  APPROXIMATIONS  OF  ODE'S 

R.  D.  Grigorief f* 

1.  Introduction.  The  main  part  in  proving  the  convergence  of  finite  difference 
approx  imat iona  for  boundary  value  problems  conaiata  in  establishing  suitable  stability 
inequalities.  One  common  method  in  proving  such  inequalities  is  to  start  from  an  often 
coaiparatively  siaple  a-priori-inequality  which  leads  to  the  deairad  stability  by  a 
perturbation  argument.  The  purpose  of  this  report  ia  to  giva  an  abstract  version  of  this 
procedure  in  the  framework  of  the  discrete  approximation  theory  introduced  by  Stuamel 
(10,11].  In  this  way  it  can  be  clearly  recognised  which  properties  of  the  problem  give 
rise  to  the  stability  inequalities.  As  a  further  consequence  one  obtains  fairly  siaple 
proofs  of  stability  results  used  only  very  indirectly  in  the  concrete  context  or  being 
even  not  found  there  (e.g.  (9,121). 

As  an  application  vs  eatablish  the  inverse  stability  of  compact  finite  difference 
approximations  for  linear  m-th  order  two-point  boundary  value  problems  on  nonunifora  grids 
in  the  maximum  norm,  which  makes  it  easy  to  get  the  convergence  of  the  schemes  obtained  in 
(9,12].  These  stability  inequalities  have  also  been  given  in  [4]  using  a  different  manner 
of  proof.  The  method  developed  in  this  report  opens  the  possibility  to  treat  also  more 
general  schemes  as  well  along  with  the  associated  eigenvalue  problem.  Moreover  the 
abstract  result  applies  equally  well  to  other  discretisation  methods,  e.g.  to  Galerkin  and 
quadrature  methods,  a  fact  which  is  mentioned  only  briefly  here  but  which  is  wellknown  to 
those  familiar  with  discrete  approximations. 
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2.  totitlaM.  Let  a  denumerable  sequence  A  and  no  reed  spaces  E,E  ,  i  e  A  ,  be 

o  i  o 

given.  The  spacea  are  said  to  fore  a  discrete  approximation  AtE.ItE^ )  if  a  linear  nap 
Lie  is  defined  on  a  linear  subset  of  all  sequence*  E]ti  e  Aq,  with  range  equal  to 

B  and  the  property 

Urn  (u,  ,l  «  A  )  ■  0  <->  lu  I  ♦  0  (ISA) 

l  O  l  o 

(the  sane  aynbol  is  usad  hare  to  denote  all  the  various  norma  on  B,  a  sequence 

ut.«  e  Aq,  lying  in  the  daauiin  of  Lie  ia  said  to  be  discretely  convergent  and  we  write 

u  ,♦  u  (i  e  A  > 

l  O 

if  u  is  its  image  under  the  nap  Lie.  By  A^,  A^  we  denote  final  sections  of 

A  ,  by  A,  A*  subsequences  of  A  ,  not  necessarily  the  sane  at  different  ocurrences. 
o  o 

The  convergence  of  a  subsequence  u  ,  i  e  A,  to  u  e  B  is  defined  in  the  obvious  Banners 

u  ♦  u  (l  e  A)s  <«>  v  ,  t  e  A  i  v  ♦  u  (i  e  A  ),  lu  -v  I  ♦  0  (i  e  A). 

I  1  O  1  Oil 

By  B  we  denote  the  closed  ball  of  radius  1  in  E  .  Let  G  ,  i  e  A  ,  be  subsets  of 
l  l  l  o 

E  .  Then  we  introduce  the  limit  set 

Lie  sup  G  s  “  (u  e  B  |  AcA  ,  u,  9  G  .  l  e  A,  u  ♦  u  (l  e  A)} . 

‘  Oil  l 

The  sequence  G(,  i  e  A^,  is  said  to  be  (locally)  discretely  compact,  if  each  (bounded) 

sequence  u(  e  G^, l  e  AcAq,  contains  a  convergent  subsequence. 

let  also  A(F,np ^ )  be  a  discrete  approxination  of  noraed  spaces.  A  sequence  of  linear 

nappinga  Lt  :  E  ^  ♦  F  ^ ,  i  e  Aq,  is  said  to  be  discretely  convergent  to  a  linear  upping 

L  >  B  ♦  F,  in  symbols  1(  ♦  1  (l  f  Aq),  if 

u  ♦  u  (i  e  A  )  «>  L  u  ♦  Lu  (»  e  A  ). 

I  oil  o 

The  sequence  L^,  t  e  A^,  is  said  to  be  consistent  with  L  if 

v  u  e  b,  u,  8  i,i  i  e  A  i  u,  ♦  u,  l,u.  ♦  Lu  (i  e  A  ). 

II  O  l  l  l  o 

The  sequence  L(,i  e  Aq,  le  said  to  be  discretely  compact  if  the  sequence 
R(L{),  i  e  Aq,  of  ranges  of  L(  is  discretely  compact.  Sometimes  we  use  the  notation 
B(,  t  e  (0),  to  Indicate  the  space  E,  and  similarly  for  F(,  L^  etc. 

3.  The  perturbation  theorem.  Let  A(B,nE(),  A(F,nr(>,  A(G,J>,  Itet<J>>,  j  -  1,2,  be 
discrete  approximations  of  noraed  spaces.  Let 
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b«  linear  mappings.  We  wish  to  conclude  that  the  stability  inequality 

<1>  ,2,utl  <  ,1«<2>uiI  *  1  u»  e  »  «  A? 

holds  if  the  inequality 

(2>  lutl  <  lL1(1)u1l  +  llll«1»u1l.  ut  e  Bj,  i  «  A1 

ia  known  to  hold.  Here  Tj«  Tj  denote  positive  constants  independent  of  u(  and  t. 
For  brevity  we  set 


Kt  «  -  Lt(2),  »  e  AoU{0)  . 

(3)  tat  (2)  and  the  following  conditions  be  fulfilled! 

<i>  lhe  sequence  t  €  A^,  is  discretely  compact  and  convergent  to  K 

(11)  The  sequence  M^Nb  )  is  locally  discretely  compact 

(lii)  The  sequence  is  consistent  with  (t*1*,  M(1)) 

(iv)  M  <2)  ♦  M(2)(t  e  A  ) 
l  o 

(V)  (w.g)  e  Lin  sup  «1(Bi),«|1)|  n  m  lop  a^11,  k^’hb^  ■> 

Hues  I  L(1)u  -  W,  Bl*1  *u  -  9 
(vi)  L(2>U  ■  0,  M*2*U  "  0  «>  U  "  0. 

Wien  inequality  (1)  holds. 

Froof .  Suppose  (1)  not  to  hold.  Wien  there  exist  a  subsequence  A  C  A  and  elements 

o 

u(  6  I  ,  l  (A,  such  that 

(4)  lUjI  -  1,  x.l(2)u,  ♦  0,  Mt<2>u,  ♦  0  (»  e  A). 

Because  of  the  assiasptlon  (i),  we  choose  A  such  that  also 

K,ut  ♦  w  <«  e  A) 

for  soaw  w  e  F.  He  now  distinguish  two  cases. 

First  caset  The  sequence  H^'^u  ,  i  (A,  is  bounded.  Wien,  since  (ii)  is  assumed, 
without  loss  of  generality 


for  some 
(5) 


•i1  * 1  *uf  ♦  g  (i  e  A) 


9  e  g 


(1) 


With  (w,9>  given  in  this  way  we  can  find  a  u  with 
9* 


r.(,,u 


„<i) 

W,  M  U 
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Now  chooso  a  consistency  sequence  y%  ,  i  e  A,  i.e. 


Replacing  Uj  in  (2)  by  ut  -  y  ,  ona  obtain* 


Yi,ui“yil  <  ,l»<2>u«  +  Ll<1)“l11  ♦  ,M,<1>(u,"*l)  '  ♦  0  <>  e  A) 


which  show*  the  convergence  u  +  u  (i  e  A).  Proa  this,  (4)  and  (i),  (ir) 


kjui  ♦  ru  -  w,  ♦  *<2>u  -  0  (i  e  A> 


and  so,  because  of  (5),  we  conclude  L*2*u  -  0,  which  gives  u  ■  0.  But  a*  a  rnn sequence 


of  (4)  it  ie  seen  u  t  0  and  we  have  reached  a  contradiction. 


ID- 


Second  case:  The  sequence  M(  u} (  l  (  A,  is  unbounded.  By  renorning  u  and  if 
necessary  passing  to  a  subsequence  skieh  we  continue  to  denote  by  A  we  obtain  a  sequence 


ut  e  E(  such  that 


(6) 


lutl  ♦  0,  Li(2>ul  ♦  0,  M1<2>uj  ♦  0,  mt(,>u  I  >1  Hi  A). 


now  proceed  in  the  saae  way  aa  in  the  first  case,  this  tias  knowing  g  f  0.  Since  g 


.CD 


H  'u  and  u  -  0  this  establishes  the  desired  contradiction. 


4.  Special  cases.  The  following  specialisation  of  theorsa  (3)  is  taylored  to  the 
treataent  of  an  ordinary  differential  operator  under  two  different  sets  of  boundary 
conditions  (see  [1,2, 4,7] ).  By  N(L)  we  denote  the  kernel  of  the  napping  L. 

(7)  lot  linear  napping* 

Ll  ;  B,  ♦  Pt,  «  Bt  ♦  S,  »  e  Ao  (0),  j  -  1,2 

be  given  such  that  with  a  constant  T1  >  0 

(8)  Tr^UjI  <  ILjU^I  ♦  lHl<1)u|l,  ut  tl|(  it  A1 
holds  and  let  the  following  conditions  be  satisfied: 

(i)  die  N(Ii)  -  din  G  <  « 

Cii)  The  restrictions  H^*|H(L),  j  -  1,2,  are  injective 

(iii>  The  eeguence  »  e  Ao,  is  consistent  with  (L.M*1*) 

ft  V  U  (2)  y(2)  .  |  . 


Proof.  We  taka  -  L/2>»  ®ls 


_  <d  _  _  (2)  „d>  _(2> 

<5 1  “  G  »  G 


G,  X  -  0,  X  -  0  In 


(3)  with  tha  discrete  convergen ee  in  A(G,0G  )  to  ba  tha  convarganea  in  G.  Hi  an  (3  >  C  ±)  — 


(iv),(vi)  hold.  Bat  alao  (3)  (v)  la  satisfied  ainea  w  •  0  and  M 
N(L)  dua  to  tha  assumptions  (7)(i),(ii). 


(1) 


la  aurjactiva  on 


Tha  naxt  specialisation  of  (3)  ia  motivated  by  tha  wish  to  atart  from  an  a-priori 


(1) 


inaquality  for  tha  oparator  L(  ,  and  ao  obtain  a  atability  inequality  for  tha  oparator 


L(*2*  diffaring  fro*  L(  * 1  *  by  ‘lower  order  tana*. 


(12)  hat  tha  following  condltlona  ba  satisfied! 


(i)  KjC  G^1*,  i  6  Aq  {0),  algebraically  and  topologically,  and  tha  aoouanca  of 


(1) 


natural  lnbeddlnga  J(  i  +  S(  ,  i  e  la  diacratalv  compact  and  convergent  to  tha 


natural  imbedding  J  i  t  ♦  6 


(1) 


(ii)  Tha  sequence  i  e  A^,  la  diacratalv  contact  and  convergent  to  K 


(1) 


(ill)  Tha  aoquenco  ,  »  «  A^,  la  conalatant  with  L 


(1) 


(iv)  Jjuj  ♦  g,  ♦  w  (i  e  A)  ->  u  e  E  i  L(t,» 


w,  Ju 


(▼> 


M  (2)  +  „(2)(l  g  A 
X  o 


(vl)  l/2*u  *  0,  M*2)u  *  0  «>  u 


0. 


■*.*.  ft*  «7£5i°ri  inaquality 
(1) 


(13)  Y^u,'  <  IL,  +  'J,®,'.  *>,  e  *,«  l  e  A,, 

holda  with 


la  T1  >  0,  than  alao  (1)  holda  with  ao*a  conatant  y2  >  0. 


Condltlona  (3X1).  (iv)-  (vi)  directly  oorraapond  to  (12X11),  (iv)-(vi).  sinca 


(1) 


J,  and  J  ,  i  e  A  ,  la  dlacretaly  contact,  (3)(ii)  ia  aatiafiad.  The 
i  I  o 


conalatancy  (3X111)  follows  fro*  (12)(iii)  and  the  convergence  J  ♦  J  assumed  in 


(12X1). 


It  ahould  be  reauurked  that  tha  condition  Ju  “  g  in  (12)(lv)  ia  only  notational  since  it 
only  diatinguiahaa  tha  element  u  C  E  and  ita  inbedding  into  G* 1  * . 

Tha  laat  apacialication  of  theora*  (3)  we  give  ia  connected  with  the  atudy  of 
dlacratiaatlona  of  aquations  of  the  second  kind  with  a  compact  operator. 
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(i)  Y,  Itt^l  <  I  itu(i,  U(  e  »  e  A,  with  BB  Y,  >  o 
(li)  L  i  I  ♦  r  is  surjective 

(ill)  The  — quence  K|#  »  e  A^,  la  dlscretslv  compact  and  convergent  to  It 
(iw)  (L-K)u  “  0  “>  u  “  0 

(v)  The  sequence  LJ#  »  e  Aq,  la  consists nt  with  L. 


Hwn.  tor 


Ya  >  0, 


Yj  iuti  <  I  cl1-«i)u1»,  u,  e  *  ,  i  e  a2  . 


Proof.  We  apply  theorem  (3)  with  L, f1)  «  “  L  ,  M,  (1 *  -  H. <2)  t  -  0,  i  e  A  {0}.  It  la 

l  ill  o 

aaay  to  chack  that  all  assumptions  of  (3)  are  satisfied. 

S.  An  application.  In  this  section  we  wish  to  show  how  tha  convergence  thaoraa  contained 
In  [9,12]  can  be  derived  frca  tha  results  of  this  note.  Incidentally,  we  establish  a 
stability  inequality  which  has  been  looked  for  in  [12,  p.  743]. 

In  [9,12]  compact,  implicit  difference  schemas  have  been  derived  for  a  ainqla  n-th  order 
differential  equation 

(16)  Lu(t)  i  -  u<B)(t)  ♦  l  a  (t)uU)(t),  t  e  [A,B] , 

l<m 

with  boundary  conditions  m 

(17)  M(q,u  «  -  U(q)[u]  +  vlq,[u]  «  -  [Vq)u(1>(A>  +  l  b(q,ull>(B)  -  C(q) 

i-0  1  i-0  1 

for  q  -  0,  . ..,  m  -  1.  These  schemas  are  of  the  form  j 

(IB)  VVV  ,  -jo  -  v(V  1 

(19)  ,  -  u‘q,(V  ♦  v‘q)[uh3  -  C(q)  ♦  c£>[fl  4  cj«‘  If] 

where  k  ■  0,  . . . ,  n-m,  q  ■  0,  . . . ,  a  -  1  and 


V  ,  .  m-1 

S'  H  tx\  *  n  T 


V  (V  •  -J0  V*  DiVA>  *  "  j0  *V*» 


(i-.u-D+ 


cXtfi ■ - j 


l  s!q)  f(t  ♦  hi.  > 

j.i  o  j ,m 


and  v‘q>.  c<q>  defined  analogously.  Here  0.  denotes  the  forward  difference  quotient 

n  V  all  • 


belonging  bo  thn  underlying  nonunifora  grid 

»h  .  -  {  bj,  j  -  °,  1 . "  |  to  -  A,  tn  -  B,  tJ+,  -tj  +  hj,  j-0, . .  -  1} 

and  H  the  distance  _,-t0  while  (j)+  means  the  positive  part  of  j.  The  meaning  of 

the  points  t  . ,  C  la  described  in  the  papers  cited. 

*«j  j#>* 

one  of  the  main  results  in  (9,12]  is  that  the  difference  operators  can  be  chosen  in  such  a 
way  that  the  reaulting  approximations  are  exact  for  polynomials  up  to  a  certain  degree  and 


the  coefficient  of  Ih»  uh»  vh> 


are  bounded  for  h  ♦  0 


where  the  normalising 


condition 

(20) 

has  been  imposed. 


J 

l  8.  .  -  1,  k  -  0,  ....  n  -  m 

j-1 

For  the  purposes  of  this  section  it  is  sufficient  to  assume 


(21)  max  {  ll^plt^)  -  I^tplt^)!,  k  -  0,  . ...  n  -  m)  ♦  0  (h*0) 


lu^lpl  -  W<q)(p]  -  C^lLp]|  ♦  0 

(22) 

Iv^tp]  -  v<q)(pl  -  (Lp]  |  ♦  0 

for  polynomials  p  of  degree  <  a  and  that 


(h  ♦  0). 


(23)  Bk,J'YS'TK'6J?J'6J?v-0‘1>'h*0- 
We  are  going  to  prove  the  following  two  stability  inequalities. 

(24)  On dor  the  conditions  (20)-(23)  there  exists  a  constant  y  >  0  such  that  for  all 


*  Z  k-0.!f?,n-i  ,Dluh<tk)l  *  *  -  ,Dq"h<to)l  +  k-O.Tff.n-m'WV1  * 

1-0  q  -  o 

(25)  Assume  (20-(23)  to  hold  and  let  (16),  (17)  be  lnlectlvs.  Then  there  exists  a 
constant  y  >  0  such  that  for  all  sufficiently  small  h  and  all  grid  functions 

k-O.^.n-i1  DiW  ‘  X  '  *  k-O.^.n-1  WV* 

i-0  q“0 
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Thar*  ar*  no  r**tr lotion*  mad*  on  th*  a*sh  ratloa  of  th*  grids  Th.  It  1*  avidant  that 
th*  convergence  r* suits  contained  in  (9,12]  ar*  an  iaaediat*  consaquanca  of  (24),  (25), 
on*  *v*n  obtain*  th*  slightly  sor*  g*n*ral  rasult  that  th*  ■  -  th  order  diff«r*nc* 
quotients  ar*  also  convergent. 

In  preparation  for  th*  proof  of  (24,  (25),  vs  r*vr it*  th*  diff*r*ne*  operator  in  th* 

fora 

(26)  Vh'V  ■  J0  \,iDiW'  k  ■  0 . .  • 

with  certain  coefficients  Inserting  successively  p(t)  “  t*,  t  -  0,  ...,  a,  into 

(21),  it  is  seen  that  (21)  is  equivalent  to  (a^  t  •  1) 

(27)  aax  {  lu^  k  “  0,  ...,  n  -  a}  ♦  0  (h*0),  i  -  0,  ...,  a. 

Wa  will  now  apply  proposition  (12)  to  prove  (25).  The  continuous  problaa  (16),  (1*  its 

into  th*  general  setting  by  taking 

E  -  c"(A,B],  r  -  C(A,B],  G(1>  -  1  [A , B)  ,  Gt2>  -  Bt* 

and  defining 

m(2)  ,  c-(.,bl  u  V  (*(o>u . e  »■. 

The  indexed  tens  ar*  defined  as  discrete  analogs.  Ms  writ*  h  instead  of  i  as 
index.  For  a  grid  function  we  introduce  th*  noras 


«  “  l  aax  {  iD^u^lt^)!,  k  •  0,  ...,  n  -  i),  *  -  a  -  1,  a. 


Then  E^,  G^1 *  are  taken  as  th*  vector  space  of  grid  function*  on  Th  equipped  with  th* 
nora  1*1  ,  I  *  1^  ^  respectively.  Fh  is  taken  to  be  th*  vector  space  of  grid  functions 
wh  defined  for  tfe,  k  -  0,  ...,  n  -  a,  noraed  by 

•whlo  i  “  aax  {|wh(tJc>l,  k  “  0,  ...»  n  -  a  ). 

(  2  )  ■ 

Finally  Gh  “  B  .  The  spaces  defined  so  far  are  easily  shown  to  fora  discrete 
approx iaatlons  (se*(6)).  Evidently,  Lj,  maps  Eh  ♦  F  , 

-(»)  -  T<2>  .  T 

\  \  *  \  * 

<2>  i  Eh  ^  ♦  (E^0^,  ...,  e  *"• 


-8' 


w«  are  now  going  to  chock  all  the  conditions  listed  in  (12).  Condition  (i)  is  proved  in 
[6,  theorem  I).  Since  K  “  0,  JC^  «  0  condition  (ii)  is  trivially  satisfied.  The 


consistency  (iii)  follows  from  (27)  and  the  fact  that,  for  each  u  e  B  and  j  *  0,  . .., 

», 

■ax  {  |D^u(ty )  -  u'^Jt^)  |,  k  ■  0,  ...,  n  -  j}  +0  (h*0). 

For  the  same  reason,  the  sequence  M^2*  is  consistent  with  M*2*.  Moreover,  due  to  the 
(2) 

structure  of  we  have 

,Mh2>uh'  <  r  ,uh'a'  “h  e  V 

with  T  independent  of  u^  and  h,  i.e.  the  stability  of  the  sequence  M^2J . 

Consistency  and  stability  together  imply  (v>.  Condition  (vi)  is  among  the  assumptions  of 
(25).  Taking  (27)  into  account  for  i  «  a  it  follows  that  |cl  |  >  a  >0  for 

KfB  O 

sufficiently  small  h  and  hence  using  the  triangle  inequality 

a  ID  u.  I  <  I  Uu.  I  +  I  y  a  ,  D.u.  I  . 
o  an  o  n  n  o  *,i  in  o 


Because  of  (27),  the  a  .  are  uniformly  bounded  and  we  obtain  for  sufficiently  small  h 

X,i 


“  'V-  <  '  Wo  +  r  'Vm-1'  S  e  S' 


o  h  m  n  h  o 

which  is  the  a-priori  inequality  (13),  It  remains  to  show  (iv).  Assume 


JhUh  *  g  (h*0)  in  A(G0>,  H  G^11) 


(28) 

(29)  I^uh  ♦  w  (h*0)  in  A(F,J!Ph). 

From  (28),  taking  (27)  into  account, 

(i) 


5  a  D.«i  *  1  <*.u  1  (h*0)  in  A(F,IIF  >. 

..•#iih,.i  n 

i<m  i<m 

Then  using  (27)  for  i  “  a  and  the  convergence  (29)  we  obtain 


r  (i) 

D  u.  ♦  w  •  }  o  u  (h*0)  in  A(F,JlF_). 

i>m  1  h 


(30)  D_ 

The  convergence  in  (28)  and  (30)  together  imply 

♦  u  (h*0)  in  A(E,!IE^) 

and  hence  from  (30)  the  desired  result  (iv),  i.e.  Lu  •  w,  follows.  Since  all  assumptions 
of  (12)  have  been  verified,  the  stability  inequality  (1)  holds,  which  is  just  (25)  in  our 


special  case 


The  proof  of  (24)  is  contained  in  th«  proof  of  (25)  if  the  boundary  conditions  are 

taken  to  be  the  initial  conditions  specified  in  (24).  of  course  for  the  initial  value 
problee  the  uniqueness  assumption  made  in  (25)  is  satisfied. 


ftcknovledgaaeat.  the  author  thanks  Professor  Carl  deBoor  for  his  valuable  suggestions 
during  the  preparation  of  this  report. 
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